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We utilize the combination of two standard trapping techniques, a magnetic trap and an optical 
trap in a Raman setup, to propose a versatile and tunable trap for cold atoms. The created 
potential provides several advantages over conventional trapping potentials. One can easily convert 
the type of the trap, e.g., from a single well to a double well trap. Atoms in different internal 
states can be trapped in different trap types, thereby enabling the realization of experiments with 
multi-component Bose-Einstein condensates. Moreover, one can achieve variations of the trapping 
potential on small length scales without the need of microstructures. We present the potential 
surfaces for different setups, demonstrate their tunability, give a semi-analytical expression for the 
potential, and propose experiments which can be realized within such a trap. 
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I. INTRODUCTION 

Trapped ultracold atomic gases serve as an ideal sys- 
tem to model many body systems and to investigate fun- 
damental questions of quantum mechanics. A paradigm 
phenomenon are the so called Josephson oscillations, 
which were demonstrated in a double well potential for 
many atoms [H-Q; more recently, even the tunneling of 
individual atoms through a barrier was demonstrated Q . 
A key ingredient for modelling new systems are novel 
ways of trapping cold atoms. A standard technique for 
trapping atoms of a single species is to use either static 
magnetic fields or a superposition of static and os- 
cillating magnetic fields leading to the so-called radio- 
frequency dressed adiabatic potentials [l0l - [T2| . Differ- 
ent species can be trapped in an optical trap making 
use of the so-called light shift [13] ■ The potential results 
from the intensity maximum (minimum) of a laser that 
is red (blue) detuned with respect to an atomic transi- 
tion frequency. By superpositions of different laser beams 
and intensity configurations one can create, e.g., optical 
lattices or double well potentials. In a more recent ap- 
proach, two lasers in a Raman configuration were used 
to trap atoms. Such a setup allows for example the cre- 
ation of optical lattices with a reduced lattice spacing 
compared to standard optical lattices [T^ |. 

A combination of optical and magnetic fields for cre- 
ating atom traps has been already described in several 
works. Even in one of the first BEC experiments a su- 
perposition of a far detuned laser and a magnetic trap 
was used to trap the atoms fl^ . In Ref. Deutsch et 
al. investigated the combination of a constant magnetic 
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field and a state dependent optical lattice that allows for 
the creation of a lattice of double well potentials. In more 
recent approaches, the superposition of radio frequency- 
fields and magnetic fields [13, [i3| or optical lattices 18 1 
were used for trapping atoms in tunable potentials. In 
this paper, we derive the potential for an atom exposed 
simultaneously to an inhomogeneous magnetic field in a 
loffe-Pritchard trap like configuration and two lasers in a 
Raman configuration. In this case, the non-trivial combi- 
nation of the magnetic and the laser fields cannot longer 
be reduced to a potential resulting from an effective mag- 
netic field. A direct consequence of this fact is that the 
potential surfaces for different hyperfine components of 
an atom do not only differ by a global factor but can 
be substantially different. Thus, it is possible to con- 
fine the different components of a multi-component BEC 
in traps of different frequencies or different types, e.g., 
one component in a double well trap and another one 
in a single well trap located at the barrier of the dou- 
ble well. Moreover, due to the availability of additional 
parameters, flexibility is gained in shaping the potentials 
compared to conventional traps. For example, one can 
smoothly convert a single well potential into a double 
well potential or drive a double well potential by vary- 
ing an offset magnetic field. Furthermore, it is possible 
to rotate the potential around one axis by changing the 
phase between the Raman lasers. 

The paper is structured as follows. In Section II, we 
derive the effective Hamiltonian. In Section III our nu- 
merical results are presented. Specifically, we provide an 
overview of the potential surfaces of the different compo- 
nents and investigate the transition of a double well to a 
single well potential as well as the rotation of a double 
well potential for one component in detail. In Section IV 
we derive a semi-analytical expression for the potential 
surfaces; in Sec. V loss mechanisms are discussed. Fi- 
nally, in Sec. VI we summarize our results and mention 
directions for possible future studies. 
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II. ANALYTICAL CONSIDERATIONS 
A. Hamiltonian and Setup 

The Hamiltonian of an (alkali) atom simultaneously 
exposed to magnetic and laser fields reads 

H = -|^Vr2 + H<=(r) + yiP(r, R) + l/^^(r, R). (f ) 

Here, R denotes the center of mass coordinate of the 
atom and r the coordinate of the valence electron rela- 
tive to the center of mass position; M is the total mass 
of the atom. H'^{r) accounts for the field-free electronic 
structure of the atom; for the scope of this work, we 
use ^''Rb as a paradigm. F'^^(r, R) and V^^^(r, R) de- 
note the contributions of the Raman lasers and the mag- 
netic field, respectively. In order to solve the coupled 
Schrodinger equation associated with Hamiltonian ([I}, 
we employ a Born-Oppenheimer separation of the center 
of mass motion and the electronic degrees of freedom. 
We are thereby led to an effective electronic Hamilto- 
nian that parametrically depends on the center of mass 
position, 

H'^iR) = H%r) + V^^{r; R) + V^^ {r; R). (2) 

Its solutions Vk(R) serve as adiabatic potential energy 
surfaces for the center of mass motion of the atom; each 
of these (trapping) potentials is then associated with a 
given internal state n of the atom. 

Regarding the magnetic field configuration, we con- 
sider the setup of a lofFe-Pritchard trap [l^ which is 
given by a two-dimensional quadrupole field in the xi , 0:2- 
plane together with a perpendicular offset (loffe-) field 
in the xa-direction; it can be parameterized as B(x) = 
GxiBi — 0x2^2 + Bje^. G denotes the magnetic gra- 
dient of the two-dimensional quadrupole field and Bj 
the constant offset field oriented along the a;3-axis. The 
quadratic term Bg oc (a;| — /9^/2)e3 that usually arises for 
a loffe-Pritchard configuration can be exactly zeroed by 
geometry, which we are considering in the following. In 
actual experimental setups, Bg provides a weak confine- 
ment also in the xa-direction. Omitting Bg, the magni- 
tude of the magnetic field at a certain position x in space 
is given by |B(x)| — + G^p^, which yields a linear 

asymptote |B(x)| Gp for large values of the coordi- 
nates {p = \/ xf + X2 ^ B/G) and a harmonic behavior 
|B(x)| w B + l^p'^ close to the origin (p < B/G). The 
magnetic field interaction within the Born-Oppenheimer 
approximation reads 

V'^ir;Il}^gpPBF-B{Il). (3) 

The Raman configuration of the excitation lasers is 
depicted in Fig. [TJ It consists of two oppositely circu- 
lar polarized lasers that are close to resonance with the 
Di transition line, i.e., being blue-detuned by A with 
respect to the transition from the bSi/2,F = 1 ground 




mp = —2 rrtF — —1 vi\.f — mp = 1 mp = 2 



FIG. 1: State linkage diagram of the unperturbed atom. 
The solid arrows denote the transitions allowed for and 
a~ light, whereas the dashed arrows indicate allowed tran- 
sitions with TT polarized light. The energy gap between the 
5-Pi/2,-F = 1 and 5Pi/2,F = 2 manifold is Ahfe = 2-n-h x 0.8 
GHz |20|]. 
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FIG. 2: (a) Setup showing the propagation direction of the 
laser beams (large arrows) and the configuration of the mag- 
netic trap. The Helmholtz coils generate the homogeneous 
lofFe field oriented along the Z direction and the loflte bars 
the quadrupole field in the X — Y plane, shown in subfigure 
(b). 

state manifold to the bPi/2,F = 1 excited state of *^Rb. 
The propagation direction of the lasers is chosen to co- 
incide with the direction of the constant loffe Field Bj. 
The overall setup is shown in Fig. [2] 

Within the dipole approximation, the potential of an 
atom exposed to the laser fields is given by 

l/^^(r;R) =-er-E(R,t), (4) 

E(R, t) denoting the electric field of the lasers. The latter 
can be expressed as 

E(R,t) = -^[E,(R,t) + E*(R,0] (5) 

i=l 

with Ej(R,i) = ejej(R)e-*('''^-'^'*+'^'(*» being the elec- 
tric field associated with the ith laser. The ampli- 
tudes £i(R) of the electric fields are spatially depen- 
dent in order to account for the focussing and shape of 
the laser beams. The factors e.^ are the unit polariza- 
tion vectors given by ei = ^(ei -f 162) (cr+ light) and 
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£2 = -^(ei — 162) (cr light), respectively. (j)i{t) take into with 
account the phases of the lasers, which additionally can 
depend on time. Expanding the atom field interaction in 
the basis \a), defined as the eigenfunctions of the field 
free atom, leads to 



F^^(r; R) = 1 ^ + c.^)|7) {a\ + h.c. (6) 



i—l CK,7 



with = e(a|E.(<)r|7) and l^I^. 
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B. Rotating Wave Approximation 

We employ the rotating wave approximation 21] in 
order to remove the time-dependence of Hamiltonian ^ 
that arises due to the laser interaction. For reasons of 
simplicity we assume that both lasers have the same fre- 
quency uji = uj and the same profile ei(R) = e(R). Since 
the magnetic field interaction term is block diagonal, i.e., 
does not mix states with different total angular momenta, 
y^^(r;R) is not affected by the transformation into the 
rotated frame. The transformed Hamiltonian reads 



RWA 



= V 



IP 



RWA 



AF 
RWA; 



(7) 



where 



^RWA ^Y.^-\a){a\+J2{Ei- MlOai (8) 

a I 

^RWA = i^5]e(R)e-(''W'(*»(?|e.r|a)|0(«| +h.c. 

(9) 



i— 1 a,/ 



(for clarity, we omit the arguments r and R of the po- 
tentials in the following). Here, a labels the different 
states of the ground state manifold and I labels the ex- 
cited states. Since the lasers are close to resonance to the 
Di transition line, we can restrict our basis to states close 
to the ground state and the first excited state. Using as 
a basis all hyperfine states of the 55* and 5P manifolds 
leads to an effective 32 x 32 matrix that will be diagonal- 
ized. 



C. Van Vleck Perturbation Theory 

In the last subsection we derived the Hamiltonian ex- 
panded in the eigenfunctions of the unperturbed atom. 
One can use van Vleck perturbation theory j22j to adi- 
abatically eliminate the excited 5P levels that serve as 
intermediate states for the Raman transitions. In this 
manner, the Hamiltonian can be reduced to an operator 
acting only on the ground state manifold (i.e., all states 
|a) = 1551/2, F = with mp G {0,±1}), 



Hyy = V 



IP 



J2Ec.\a){a\+J2^Po.\l3){c 



(10) 



1 



'a — El Ej3 — El 



(11) 



being the effective interaction within the ground state 
manifold. Here, the index I labels the excited states, 
which have been eliminated. For a detailed derivation of 
Eq. pT|) . we refer the reader to the appendix of Ref. [23| . 
Employing 

-\jZY.<^)^"'^'''''^''^'^Hl\e.v\a). (12) 

i—l a,l 

for the block off-diagonal matrix elements of V[^^ yields 
the effective interaction W within the ground state man- 
ifold whose matrix representation correspondingly reads 

1 2 



x( 



i.i' — l I 

1 



Ea - Ei+ huj Ep - Ei + hw 



(13) 



If one restricts the sum over the intermediate states 
to the 5Pi/2 hyperfine sublevels, which represent a good 
approximation, one obtains the compact form 

1(5^1/21 |er||5Pi/2)P / ^ C\ 
^ ^ \\ — 1/211 II — iJ2!}_ B (14) 

72ceoA j^j 

for the atom laser interaction, with 



A 



B = 2\l 



C 




/(R), 



/(R), 



e^^'^(*)/(R) 



(15) 
(16) 
(17) 



A = -5551/2,^^=1 ~ -^5Pi/2,F=i - being the detun- 
ing of the transition lasers and A$(i) — (f)i{t) — 02 (i) 
their phase difference. /(R) = ceo|e(R)P/2 denotes 
the intensity of the laser, being the dielectric con- 
stant and c the speed of light. The reduced matrix el- 
ement in Eq. (|14p can be deduced from the measured 
lifetime of the excited state which yields in our case 
(55i/2||er||5Pi/2) = 2.9919eao [H, Hi. The individual 
contributions of the matrix (|14p can be interpreted as 
follows. The diagonal elements stem from the light shift 
potential of the lasers, i.e., the off-resonant coupling of a 
mp component of the ground state to an excited state. 
The off-diagonal elements arise due to the coupling of 
the mp — I {mp = — 1) component via an intermediate 
(excited) state to the mp = —1 {mp = 1) component. 
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These off-diagonal matrix elements are not present in ra- 
dio frequency traps. The specific form of the matrix oc- 
curs since the laser light is circularly polarized. Other 
polarizations would lead to a coupling of different states, 
i.e., different off-diagonal entries of the matrix W. Fur- 
thermore, we should note at this point that the matrix is 
expanded in the field-free basis of the atom. In this ba- 
sis, the contribution of the magnetic field interaction V^^ 
for the ground state manifold is represented by the spin 
matrices for total spin F = 1, giving rise to off-diagonal 
matrix elements as well. As a result, the combined action 
of the magnetic and laser fields leads to the mutual cou- 
pling of all magnetic sublevels of the 55*1/2, = 1 man- 
ifold, represented by a fully occupied matrix. In Section 
IIVI we tackle this issue by performing a principal axis 
transformation that diagonalizes the magnetic field in- 
teraction and thus provides us with a more suitable basis 
for the interpretation of the underlying physics. 



III. NUMERICAL RESULTS 

We restrict our investigations to both lasers having a 
Gaussian profile with a width a in the X, Y plane and 
assume that they have a constant intensity in the prop- 
agation direction Z , 



J(R) = lo exp(- 



(18) 



Since the magnetic field interaction term is independent 
of Z as well, we find a total potential that is constant 
in the Z-direction. In order to provide a confinement in 
the Z-direction one can, e.g., utilize an additional laser 
or make use of the defocusing of the laser beams. If 
not stated otherwise, we fix the detuning of the lasers 
to A = — Aiifs/2, i.e., right in the middle between the 
5Pi/2,F = 1 and the 5Pi/2,F — 2 excited states as de- 
picted in Fig.[Tl and the intensity to /o = 10 W/m^. 



A. Overview Over All Components 

Diagonalizing the Hamiltonian matrix (ITO|) leads to the 
adiabatic potential surfaces (R) for the center of mass 
motion as a function of the center of mass coordinate R. 
As expected from the magnetic field interaction, we find 
one trapped and one anti-trapped component, according 
to the quantum numbers mi? = +1 and mp = —1, re- 
spectively. Interestingly, the mp — component - that 
is untrapped in a pure loffe-Pritchard field - now shows 
an attractive potential with a double- well structure along 
the X-axis. Note that the quantum number mi? is only 
valid in a rotated frame of reference that we are going 
to introduce in Sec. |lVl Nevertheless, we continue to use 
mi? as a label for the different states in the laboratory 
frame in order to avoid confusion. Since the different po- 
tential surfaces are separated, one can uniquely assign to 
each state one particular surface. 



3. 




Position [fim] 

FIG. 3: Potential curves for Y = (solid line) and for X — 
(dashed line). Two components show an attractive potential 
in two dimensions whereas one component is exposed to a 
repulsive potential. The individual potential surfaces are well 
separated. 




-10 

Position [fjva] 

FIG. 4: Difference of the potential differences for adjacent 
components for X = (solid line) and Y = (dashed line). 
The deviation between the transition frequencies of adjacent 
potential surfaces is non-zero and depends on the center of 
mass position of the atom. 



For investigating the energy spectrum further, let 
us define the radio transition frequencies (5_i(R) = 
14.^=_i(R) - Kn^=o(R) and ^i(R) = 14.r=o(R) - 
Knp=i(R) between the mp = and the mi? = =Fl states, 
respectively. Figure |4] provides a measure for the devia- 
tion of both transition frequencies by showing the differ- 
ence (5_i(R) - (5i(R) for X = (solid hue) and F = 
(dashed line). The fact that the transition frequencies 
do not coincide, i.e., the difference being non-zero, can 
be used to mutually couple two components without cou- 
pling to the third component. This allows for example a 
transfer of atoms from the mp = —1 component to the 
mp = component without coupling to the untrapped 
mp = 1 component. Moreover, the radio transition fre- 
quencies depend on the center of mass position R and 
thus on the absolute value of the potentials. Therefore, 
one may couple energy-selectively one component to an- 
other one, i.e., couple atoms of the mi? = or mi? = —1 
component at a certain position to the untrapped mp = 1 
component. Such a scheme can be used to evaporatively 
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cool the mp = or the mp = —1 component [2i 



B. mp = Q Component 



laser intensity. Hence, one can increase the height of the 
barrier by increasing the laser intensity. The position of 
the local minimum can be varied by changing the ratio 
B/G or the width of the laser cr. 



We start our detailed investigations of the individual 
components with the potential surface for the mp ~ 
component. Note that this component is not confined 
in a pure magnetic trap. Figure [5] shows the contour 
plot of its trapping potential for A$ = 0, cr = 10 /im, 
G ^ O.lG/iim, and (a) Bi = 10 G, (b) = 1 G, (c) 
Bi 0.5 G, and (d) Bj = 0.1 G. The shape of the po- 
tential correspondingly changes from (a) close to being 
rotationally symmetric to (b) a cigar shaped potential 
in the X-direction, and finally to (c,d) a double well in 
the X-direction. Thus, one can change the shape of the 
potential from a single- to a double well potential by 
changing the magnitude of the loffe field. Alternatively, 
one may also drive the double or single well potential by 
modulating the magnitude of the loffe field. The param- 
eters of the double well trap can be tuned in different 
ways. The height of the barrier can be tuned by chang- 
ing the ratio ^ — G/B. This is shown in detail in Fig. 
[9] in Sec. IIVI The position of the minima can be con- 
trolled by changing the width of the lasers. This leads 
at the same time to a change of the height of the barrier 
and consequently to a change in the number of trapped 
states within each well. For a more detailed discussion of 
the properties of the double well potential, we refer the 
reader to Sec. IIVI where a semi- analytical expression for 
the potential surface is derived. By displacing the center 
of the laser beams in the X-direction with respect to the 
loffe Pritchard trap, one can create in addition a tilted 
double well potential. 

Figure [S] shows the effects of a phase difference be- 
tween the two excitation lasers on the potential surface. 
A phase difference A$ leads to a rotation of the whole 
potential surface about the Z-axis by A$/2. For a zero 
phase difference one can add the electric fields of the 
lasers, resulting in an effective electric field which is po- 
larized linearly along the X-axis, whereas a phase differ- 
ence of TT leads to an effective electric field which is polar- 
ized linearly along the y-axis. In general, a phase differ- 
ence of A$ leads to a rotation of the polarization vector 
of the total electric field by A$/2. The sensitivity of the 
alignment of the double well potential on the orientation 
of the polarization vector seems at first glance surprising, 
given the azimuthal symmetry of a pure loffe-Pritchard 
trap. It is rooted in the spatially varying quantization 
axis of the loffe-Pritchard trap. This issue is analyzed 
further in Sec. IIVI to which we refer the reader at this 
point. 

Figure [7] shows the potential surface for the same setup 
as in Fig. [5] but for the case of a single cr^ -polarized 
laser instead of a pair of /a~ polarized lasers. The 
potential is rotationally symmetric with a local maximum 
at the origin. Therefore it can be used as a ring shaped 
trap. The absolute value of the potential scales with the 



C. rap = — 1 component 

The potential of the mp ~ — \ component is dominated 
by the contribution of the magnetic field and therefore re- 
sembles the attractive potential of a pure loffe-Pritchard 
trap, namely a single well with a minimum at the ori- 
gin. Without lasers, the confinement in the X- and In- 
directions is equal, leading to an isotropic potential. The 
contributions of the lasers break this symmetry, giving 
rise to a slightly ellipsoidal potential. The value of the 
eccentricity depends on the intensity of the lasers. A 
phase difference between the lasers leads to a rotation of 
the (anisotropic) potential surface about the Z-axis. 

IV. SEMI-ANALYTICAL INTERPRETATION 
OF THE POTENTIAL 

A. Principal Axis Transformation 

In the previous sections we investigated the system in 
the laboratory frame of reference where the quantization 
axis for the atom is determined by the direction of the 
constant loffe field. However, because of the inhomogene- 
ity of the magnetic field, a more adequate description of 
our system is to define the quantization axis along the 
local magnetic field vector B(R). In this chapter, we 
tackle this issue by introducing the spatially dependent 
unitary transformation 

Ur = exp(-iai^x) exp(-i/3Fj^) (19) 

that rotates the local magnetic field vector into the z- 
direction of the laboratory frame of reference with the 
total spin vector F = L + S + I consisting of the sum of 
the electronic orbital angular momentum vector L, the 
electronic spin vector S, and the nuclear spin vector I. 
The corresponding rotation angles are defined by sin a = 
^GY/y/B-^ + G2(X2 + y2)^ sin/3 = GX/VW+WX^, 
cos a = V-B2 + G2X2/^B2 + G'^{X'^ + Y"^), and cos /? = 
Bj^jB'^ + G2X2. This rotation diagonalizes the mag- 
netic field contribution in Hamiltonian ([5]), giving rise 
to 

UrV'^Ul = .gFAisi^.|B(R)|. (20) 

Note that in the rotated frame of reference without lasers 
mp remains a good quantum number even in the presence 
of the inhomogeneous loffe-Pritchard field. In absence of 
the Raman lasers, the trapping potentials correspond- 
ingly read K = 5fMbWf |B(R)|. 

In order to solve the time-dependent Schrodinger equa- 
tion associated with Hamiltonian ([2]), the Hamiltonian 
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FIG. 5: Contour plot of the potential surface for the mp = component for A$ = 0, ct = 10 /im, G — 0.1G//im, and (a) 
Bi = 10 G, (b) Bi — 1 G, (c) 5/ = 0.5 G, and (d) Bi = 0.1 G. The grey coded values of the potential are given in nK. 
By modulating the magnitude of the loffe field one can transform the potential smoothly from a single well to a double well 
potential. 
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FIG. 6: Contour plot of the potential surface for the mp = 
component for G = 0.1 G/fim, a = 10 fim, Bj = 0.1 G, and 
for (a) A* = 0, (b) A$ = tt/2, (c) A$ = tt, and (d) A$ = 
37r/2. A phase difference of Acl? between the Raman lasers 
leads to rotation of A$/2 of the potential surface around the 
^-axis. 



for the atom in the loffe-Pritchard trap and the laser 
interaction must be expressed in the same frame of ref- 
erence. Hence, the unitary transformation Ur must be 
applied to V-^^ as well. We find 
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FIG. 7: Potential surface for the case of a single a polarized 
laser for a = 10^ m, G = 0.1 G/fim, and Bi =0.1 G. The 
potential is ring shaped. 



picted in Fig. [T] become 



e± ■ UrvU^ — — p [x cos l3 + y sin a sin /3 
v2 

— zcosasin/3 ± i(j/cosa -|- zsina)] . 

(22) 

Equation (j22p can be rewritten in terms of the polariza- 
tion vectors e± and eo defined in the rotated frame of 
reference, showing that in a lofFe-Pritchard trap contri- 
butions of all polarizations emerge away from the trap 
center [l^, [23] . This changes drastically the simple tran- 
sition scheme caused by the and a~ light as depicted 
in Fig. [1] leading to a spatially dependent coupling be- 
tween the involved ground- and excited states. 



/ X cos /3 -f y sin a sin /3 - z cos a sin /? \ ^- Effective potential for the mp = component 

UyvUl = y COS a + z sin a 

y a;sin/3 - y sin a cos /3 -I- z cos a cos /3 j For our setup as depicted in Fig. [T] and zero relative 

(21) phase A<f> = 0, the operator describing the interaction 

Therefore, the and laser transitions that are de- between the laser field and the atom reads in the rotated 
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frame of reference 

diX,Y) = {e++e-)-UrrU^, 

sin a sin p — \/2z cos a sin p. 

(23) 

Note that the dependence of 0{X,Y) on the center of 
mass coordinates X and Y stems from the dependence 
of the rotation angles a and /? on these coordinates. One 
finds the foUowing hmits. At the origin, the transforma- 
tion is given by unity, providing 0(0,0) = -\/2a;, i.e., tt- 
polarized Ught in the x-direction. For X = 0, 0{X,Y) is 
invariant under Ur which results in 0(0, F) = 0(0,0) — 
^J2x. For y = and X ^ oo the gradient field dom- 
inates and the operator corresponds to the operator of 
TT-polarized light in the z-direction O((X),0) = \/2z. A 
non-zero relative phase between the lasers leads to a dif- 
ferent polarization of the total electric field in the labo- 
ratory frame, giving rise to a rotation of the operators 
in the rotated frame. This explains the rotation of the 
potential as seen in Fig. [6] Since it is straightforward to 
generalize our results to non-zero relative phases, we will 
restrict our analytical considerations to a zero relative 
phase in the following. 

Within van Vleck perturbation theory as introduced 
in Sec. [H the effective interaction in the rotated frame of 
reference becomes 

W;3„=ge(R)' ^(/3|e.C/,rC/t|/)(/|e,C/,rf/t|a) 

^ ^ Ec - El + hoj ^ Efs - El + hJ' 

cf. Eq. [T31 Note, that we assumed once more that the 
shapes and frequencies of both lasers are identical. We 
are interested in the regime of large enough magnetic 
fields (adjustable by the homogeneous loffe field com- 
ponent Bi) where the Zeeman splitting overcomes the 
light shifts of the laser fields, i.e., IW/JqI ^ g_F^B|B(R)|. 
Hence, we can approximate the fully occupied effective 
interaction matrix V^^ -I- W by omitting the off-diagonal 
matrix elements yV^Q,,a ^ /3, that couple the Zeeman- 
splitted mp components. This procedure leaves us with 
the diagonal matrix elements of the effective interaction. 



U = ie(R)2^ [cos2/3|(a|. 



+ sin2asin2^|(a|y|0P 

+ cos^ a sm^ f3\{a\z\l)\^ /{E^ - Ei + Hoj). 

(25) 

Employing |(a|2:|Z)p — \{a\y\l)\'^ and performing the sum 
over all intermediate states \l) eventually yields for the 
mp = component the effective potential 



X' 



X^ + Y^ 



[Voo - V^), 



(26) 



where 



36ceo 



ihfs 



|(5S'i/2||er||5Fi/2) 



/(R) 1(5^1/21 |er||5Pi/2) 
9ceo A + Ahfs 



(27) 
(28) 



are the light shifts at the origin and in the limit Y = 
0,X — >• oo, respectively. ^ = B/G is a length scale 
characterizing the particular configuration of the loffe- 
Pritchard trap. Numerical comparison of the effective 
potential with the corresponding eigenvalue of the 
full problem (O shows a excellent agreement for small 
laser intensities. We observe a maximal relative devia- 
tion of less then l%o for / = 10 W/m^. For larger in- 
tensities, the agreement gets worse since the off-diagonal 
matrix elements Wpa increase in magnitude. However, 
even for / = 100 W/m^ the deviation is less than 5%o. 



C. Discussion of the Effective Potential 

The analytical prediction of the effective potential Kff 
for the mi? = component allows us to deduce its basic 
properties by a simple analysis of Eq. ((26)) . The effects 
of the magnetic field creating the double well potential 
described in the previous section, are included in the sec- 
ond term in Eq. (1^51) . In order to obtain the double well 
structure this term needs to be negative. This is the case 
for laser light which is red detuned with respect to the 
F — 2 state and blue detuned with respect to the F = 1 
state, i.e, > A > —Ahfs. An interesting case occurs for 
A — — Ahfs/4. Then the offset potential Vq vanishes and 
the height of the barrier is equal to the maximal depth of 
the potential. For > A > —Ahfs/4 the barrier-height 
is larger than the depth of the wells with respect to the 
continuum and for Ahfs < A < —Ahfs/4 the height of 
the barrier is smaller than the depth of the wells with 
respect to the continuum. 

In the following, we restrict our investigations again to 
the case A = — Ahfs/2. Since the double well only oc- 
curs for > A > —Ahfs, this choice leads to a maximal 
detuning with respect to the F = 1 and F = 2 states. 
As expected from the numerical solutions provided in 
Sec. (|IIII) . VcS shows a double well structure that is cen- 
tered at the origin. The positions {Xo,Yo) of the two 
local minima are given by 



Xo^± 




(29) 



and Fo = 0. It is obvious from Eq. p9|) that the double 
well only exists if the discriminant is positive, giving rise 
to the condition ^ < cr. For £, > a one finds a single well 
potential, whereas for decreasing ^ < cr the double well 
starts to build up. Starting in the limit ^ — >■ from a 
double well with a barrier with finite height but width 
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2 4 6 8 f[^m] 2 4 6 8 a [/mi] 



a more narrow laser entails a more shallow double well. 
Since the range of ^ is determined by the condition ^ < 
cr, a more narrow laser necessitates furthermore tighter 
magnetic traps, i.e., smaller ^. Note that the positions of 
the minima do not depend on the intensity of the lasers. 
Thus, by increasing the laser intensities one can increase 
the height of the barrier without changing the position 
of the minima. In this way, the number of trapped states 
in each well can be controlled. 



FIG. 8: (a) Contour plot of the effective potential for Y = 
for fixed cr = 10 /im as a function of f . With increasing ^ the 
barrier gets lower and broader, (b) Same potential for fixed 
^ = 1 and as a function of a. For a > ^, the shape of 
the barrier close to the origin is almost conserved. However, 
the position of the minima increases with increasing cr. In 
both subfigures, the dashed line indicates the positions of the 
minima. 
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FIG. 9: Dependence of the height of the barrier on ^ for 
different a. For increasing ^ the height of the barrier decreases 
monotonically. 



going to zero, the distance between the minima increases 
with ^ up to a local maximum AX^'^^ = 2{y/2 — l)a at 

Cmax = V 3\/2 — 4(7 and decreases for even larger ^ up to 
^cr — c where the minima vanish, thus transforming the 
energy surface to a single well potential. This behavior 
is illustrated in Fig. [SJa) where a contour plot of the 
effective potential along the X-axis as function of ^ is 
shown for a fixed width u = 10 /xm of the lasers; the 
dashed line indicates the positions of the minima. Figure 
[HJb) shows a similar contour plot of the effective potential 
along the X-axis, but now as a function of a for fixed 
^ = 1 /im. In this case, the shape of the barrier close to 
the origin is approximately conserved but the position of 
the minima increases with increasing a. 

The height of the barrier is given by AV = I^g (0, 0) — 
Vef[{Xo,0) and is directly proportional to the intensity 
of the lasers. Figure [HI shows the dependence of AV on 
the parameter ^ for different values of a and fixed laser 
intensity / = 10 W/m^. The behavior of AV for different 
a is qualitatively very similar: the height of the barrier 
decreases with increasing ^ monotonically. For fixed ^, 



D. Effective Potential for a Single Laser 

One can use the same semi-analytical procedure as 
above in order to predict the effective potential for a sin- 
gle laser. One finds 



1, 



off 



e 



i?2 



(30) 



where = X'^ + Y^. The comparison of the effective 
potential (j30p and the corresponding numerical solutions 
shows again an excellent agreement. As opposed to the 
Raman setup involving two lasers, in the case of a single 
laser the resulting trapping potential for the mp = 
component is rotationally symmetric. The ring-shaped 
minimum is located at 



1 



6 



6 



(31) 



and exists for a > Starting at a barrier with ar- 

bitrarily small width with maximal height for ^ — >■ 0, 
the barrier at the origin gets lower as f increases and 
eventually vanishes as ^ = c/a/Z. The distance of the 
minimum to the origin increases with ^ up to the lo- 
cal maximum at 



^max = Y 5^/3/2 — 6(7 with a value 
-Rmax = (V3 — \/2)cr. Then it decreases again and be- 
comes zero at ^cr = 1/ V^cr thereby transforming the po- 
tential into a single well. Note that the potential is not 
a simple superposition of a magnetic single well poten- 
tial and a repulsive potential created by the laser. The 
spatial structure is a direct consequence of the spatially 
dependent light shift potential of the laser, giving rise to 
a barrier which is smaller than the width of the laser. 



V. LOSS MECHANISMS 

In the previous sections we have shown that for appro- 
priate parameters the discussed combination of external 
fields leads to a confinement for two components of the 
ground state manifold. Let us discuss in this section pos- 
sible loss mechanism for these potentials. 
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A. Lifetime of the Intermediate State 

One loss channel results from the coupling of the 
ground manifold of states to the excited 5Pi/2 states. 
Despite the fact that the Raman lasers are detuned with 
respect to the excited state, there is a finite probabil- 
ity to excite the atom to this state due to the width of 
the state and the width of the lasers. The excited atom 
can subsequently decay spontaneously to the untrapped 
ground state. The resulting lifetime of the dressed state 
can be estimated by applying perturbation theory, lead- 
ing to Tcff — t(A/wcp)^ with T being the lifetime of the 
unperturbed excited state and Wcp the coupling matrix 
element of the ground state to the excited state. For our 
parameters (A — — Aiifs/2, Ii = h = 10 W/m^, and 
T = 27 ns) we get an effective lifetime of Toff ~ 27 ms. 
The latter can be increased by decreasing the intensities 
of the lasers. However, one has to bear in mind that this 
will reduce the depth of the trapping potentials as well. 

B. Inelastic Collisions 

Another loss mechanism occurs if more than one atom 
is loaded into the potential due to the mutual interac- 
tion of the atoms. This mechanism can be estimated 
on a mean field level, incorporating an effective coupling 
coefficient that determines the interaction between the 
atoms. In the field-free case one gets a population trans- 
fer from one component to another due to interaction 
when there is an overlap between the wave functions of 
two components. For our setup, however, one obtains 
new dressed states that are superpositions of the field- 
free states. Hence, one obtains a state changing contri- 
bution due to interaction even if only one dressed state is 
occupied. However, for the above discussed parameters 
this additional term can be neglected. 

VI. CONCLUSIONS AND PHYSICAL 
APPLICATIONS 

We investigated the trapping potentials for ^'''Rb 
ground state atoms simultaneously exposed to a magnetic 
trap in a loffe-Pritchard like configuration and an optical 
trap in a Raman setup. The Raman lasers were detuned 
between the two excited hyperfine states, the F — 1 

and the F = 2 state. By varying the offset field of the 
loffe-Pritchard trap, we demonstrated that the trapping 
potential of the mp = component can be tuned from a 
rotationally symmetric single well to a double well trap; 
in the intermediate regime, one finds a cigar shaped trap- 
ping potential. By applying a phase difference between 
the two Raman lasers, the resulting trapping potentials 



can be rotated about the propagation direction of the 
laser beams. A semi-analytical formula for the poten- 
tial surfaces has been derived. All relevant properties of 
the double well potential have been determined analyti- 
cally as a function of the various trap parameters. For a 
single excitation laser, the proposed scheme results in a 
ring-shaped trap for the nip = component. 

In order to exploit the unique features of the above dis- 
cussed potentials, one can think of various experiments. 
For example, one might trap the mp = component in a 
double well potential which gives for the mi? = — 1 com- 
ponent a single well potential located at the center of the 
barrier. For an asymmetric occupation of the wells one 
can then observe tunneling of atoms trapped in the dou- 
ble well potential through the atoms trapped in the single 
well potential. The oscillation frequency of the tunnel- 
ing can thus be investigated as a function of the occu- 
pation number of the second component located at the 
barrier, which is reminiscent of a single atom transistor 
[28j . We performed corresponding numerical simulations 
of the Spinor Gross Pitaevskii equation and observed in- 
deed an increase of the oscillation period with increasing 
occupation number of the second component. 

In a similar setup as mentioned above (one component 
is exposed to a a double well potential with a narrow bar- 
rier and one component to single well potential centered 
at the position of the barrier) one might trap atoms in 
the component exposed to the single well potential and 
then transfer all atoms by an rf-pulse in the "double- well" 
component that feels a sharp potential maximum at this 
point. Depending on the energy of the atoms the conden- 
sate wave function would consequently either split into 
two parts or the whole wave function would move into 
one direction, which could be used as a test of the valid- 
ity of the Gross Pitaevskii equation [2^. Moreover, the 
possibility to rotate an anisotropic single well trap (which 
arises for the mp = — 1 component) opens up a new pos- 
sibility to study superfluids under rotation [30|-|32|. The 
possibility to rotate the double well potential allows one 
to create an effective ring potential by rotating the po- 
tential fast enough so that the atoms feel a time averaged 
potential 33j . In this way the transition of a double well 
to a ring-shaped potential can be investigated. Moreover, 
one can exploit the feature that one can drive the poten- 
tial surface of the mp — component by modulating the 
loffe field st reng th to investigate non-equilibrium driven 
systems [sl. IsSj. 

As an extension to the present work, it would be in- 
teresting to extend the studies to different propagation 
directions of the laser beams with respect to the orien- 
tation of the magnetic trap. Furthermore, one can vary 
the shape of one or both laser beams by using excited 
modes. 
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